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Rossby waves drifting in the azimuthal direction are a common feature at the onset of thermal convective instability in a rapidly rotating 
spherical shell. They can also result from the destabilization of a Stewartson shear layer produced by differential rotation as expected in 
the liquid sodium experiment (DTS) working in Grenoble, France. 

A usual way to explain why Rossby waves can participate to the dynamo process goes back to[Bussc (197§|). In his picture, the flow 
geometry is a cylindrical array of parallel rolls aligned with the rotation axis. The axial flow component (the component parallel to the 
rotation axis) is (i) maximum in the middle of each roll and changes its sign from one roll to the next. It is produced by the Ekman 
pumping at the fluid containing shell boundary. The corresponding dynamo mechani sm can be ex plained in terms of an a-tensor with 
non-zero coefficients on the diagonal. It corresponds to the heuristic picture given by Bussc (1975). 

In rapidly rotating objects like the Earth's core (or in a fast rotating experiment), Rossby waves occur in the limit of small Ekman 
number (ss 10~ 15 ). In that case, the main source of the axial flow component is not the Ekman pumping but rather the "geometrical 
slope effect" due to the spherical shape of the fluid containing shell. This implies that the axial flow component is (ii) maximum at 
the borders of the rolls and not at the centers. If assumed to be stationary, such rolls would lead to zero coefficients on the diagonal of 
the a-tensor, making the dynamo probably less efficient if possible at all. Actually, the rolls are drifting as a wave, and we show that 
this drift implies non— zero coefficients on the diagonal of the a-tensor. These new coefficients arc in essence very different from the ones 
obtained in case (i) an d cannot be interpreted in terms of the heuristic picture of IBussj IU975T ). They were interpreted as higher-order 
effects in Bussc (1975). In addition we considered rolls not only drifting but also having an arbitrary radial phase shift as expected in 
real objects. 

Keywords: Dynamo effect, mean field electrodynamic, waves, Ekman number 



1 Introduction 



Rossby waves naturally result from thermal convection instabilities in a rapidly rotating shell. Different 
configurations have been studied depending on whether the fluid lies between two concentric spherical 
shells or inside a full sphere, and on the type of heating (either differential or internal). At the instability 
onset, the motion takes the form of rolls aligned with the axis of rotation and localized at the vicinity 
of a cylinder lying in the bulk of the fluid. These rolls drift usually as a wave in the prograde azimuthal 
direction. In each roll in addition to the horizontal flow, that is, the flow in a plane perpendicular to the 
axis of rotation, there is also an axial flow (component parallel to the rotation axis) due to the boundary 
conditions at the ends of the rolls. 

In rapidly rotating objects like the Earth's core, the Ekman number E, defined as the ratio of the 
viscous to the Coriolis forces, is small (E oc 1CP 15 ). Rossby wave as a linear solution of the t hermal 



conve ction problem in the asymptotic limit of small Ekman number E <C 1, was first proposed by Bussc 



confirmed numerically by Dormy et al\ ( 20041 ) (see also Zhang! 1991 . 19921 . Zhang and Jones 19931 ). These 
results assume the asymptotic limit Pr/E — >• oo where Pr is the Prandtl number defined as the ratio 



( 1970) (see also the contribution by Roberts! 1968). After several i ntermediate improvements ( Soward 
19771 . lYanollTflfli ). an exa ct solution w a s give n by I Jones et al\ (120001) . Since then, the solution h as been 
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of the visco sity to the t hermal diffusivity. Additional issues have been addressed, i n the asymptotic limit 
Pr/E — > ( Zhang) 1995 ) and in the general case < Pr/ E < oo (Zhang et aZ.ll2Q07l ) . including discussions 



about the nature of Rossby waves versus inertial waves (IBusse et al 
Thermal Rossby waves have been studied also experimentally 



200fii ). 



( Busse and CarriganI 1976 

Carrigan and Bussa 19831 . Cardin and Olson 19941 . see also the review paper by Bussd [20021 and refer- 
ences therein). Above the onset the flow becomes highly turbulent and the non-linearities may be s trong. 



Bussd 



Though depending in a comp lex way on the parameters of the problem (jGrote and Bussd 12001 
20021 . iMorin and Dormvll2004l ). it is worth noting that the per sistence of the columnar structure of Rossby 



waves has been observed both experimentally and numerically (lAubert et a/.ll2003l . lCardin and Qlsonlll994l . 
Sumita and Olson! l2000h . 

When the fluid is electrically conducti ng such Rossby waves, in combination with differential rotation, 
are expected to pro duce dynamo act i on ( Kagevama and Sato 19971 ). In addition, by processes related to a 



2D in verse cascade (Sommerialll986l . lAubert et al. 2001 ) or to the presence of a strong toroidal magnetic 
field ( Cardin and Qlsonlll995l ). the number of rolls for a very low Ekman number is expected to be much 
lower than estimated from the asymptotic theory of the onset of thermal convection. Therefore, though 
the parameters in numerical simulations or experiments are far from those of the Earth's outer core, these 
studies suggest that the existence of Rossby waves in the form of columnar structures of reasonable size 
may occur and be important in the geodynamo process. 



Ros sby w aves can also be obtained me chanically, instead of t hermal ly, as shown by iHide and Titman 
(11967m and |Bussi (jl968l ). More recently ISchaeffer and Cardinl (|2005al ) considered a fluid between two 
concentric spherical shells in fast rotation but with slightly different rotation rates. They found that the 
destabilization of the Stewartson shear layer at the tangent cylinder leads indeed to Rossby waves. They 
have shown that such Rossby waves with th e strong differential rotation present in the fluid are capable of 
dynamo action (jSchaeffer and Cardinl 120061). Thi s has enforced the interest in building a new experiment 
in liquid sodium, called DTS (jCardin et qiJlioO^ . iNataf er^1l2006l . ISchmitt et fllll2008h . 

In dynamo theory, it is known for long that differential ro tation, that is the f2-effect, can generate a 
toroidal magnet i c field from a poloidal magnetic field (see e.g. lElsasser 1956). The comprehensive picture 
given by Bussd (jl975l ) illustrates how the interaction of Rossby waves with a toroidal magnetic field 
can generate toroidal electric currents and thus a poloidal magnetic field. Within the mean-field concept 
this process can be described in terms of an a-effect. The overall dynamo mechanism is known as an 
af2-dynamo. 

We can show by simple arguments, that the efficiency of the process described by Bussd (jl975l ). or 
the magnitude of the a-effect, depends critically on the relative positions of the horizontal and axial 
components of the flow in each roll. As a first step let us ignore the drift of the rolls. As illustrated in 
figure 12k, each roll rotates around its axis with a rotation rate changing its sign from one roll to the 
adjacent ones. We distinguish between the two cases in which the axial flow is maximum either (i) within 
each roll or (ii) between two adjacent rolls (top row of figured]). A given large scale azimuthal magnetic 
field B is stretched by the fluid motion leading to a secondary magnetic field b (middle row of figure [TJ . 
Then an azimuthal electromotive force u x b is created (bottom row of figure [TJ). For case (i) all local 
electromotive forces act in the same sense and so generate a global azimuthal electric current. This can be 
interpreted as an a-effect. On the other hand, for case (ii) local electromotive forces with opposite signs 
occur, implying eventually that there is no global azimuthal electric current, that is, no a-effect. 

In a rapidly rotating shell with a rigid boundary for which the no-slip condition applies, and relying 
on the quasi-geostrophic approximation, the axial flow u z in each roll is the sum of two terms. One term 
is the Ekman pumping, scaling as E 1 ! 2 and its intensity is maximum within each roll as in case (i). The 
other term is the geometrical slope effect at the ends of the rolls, scaling as E° and proportional to the 
ra dial flow u K . It is the n of ma ximum intensity between two adjacent cells as in case (ii). It is argued 



in 



Schaeffer and Cardinl (|2005ah that in rapidly rotating spherical shells like the Earth's liquid core, the 
Ekman number is so small that the first term in u z can be neglected compared to the second one. Then 
in the light of the arguments illustrated in figure [1] (ii) the ability of such a flow to generate an a-effect 
for £ C 1 could appear questionable. Actually in such rapidly rotating systems, stationary convection 
can not occur. The rolls have to drift as a wave. Then taking this drift into account, we shall demonstrate 
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Figure 1. In two adjacent rolls, the fluid velocity u (small arrows and inwards/outwards symbols in first row) acting on a given azimuthal 
mean magnetic field B (thick arrow in top row) results in a secondary magnetic field b (middle row), and so in an electromotive force 
u X b (bottom row). The radial, azimuthal and axial directions correspond to s, ip and z. Left column: case (i), for which the axial flow 
is maximum within each roll. Right column: case (ii), for which the axial flow is maximum between each roll. In case (i) u and b are 
clearly not parallel or antiparallel, and the azimuthal component of u X b has the same sign everywhere in the considered points, that is, 
its azimuthal average is non— zero. By contrast in case (ii) u and b are presumably more or less parallel or antiparallel in the considered 
points so that uxbis small. Moreover the azimuthal component of u X b has different signs in each cell. These aspects suggest that its 
average vanishes. 



th at even in ca se (ii) an a-effect remains possible. A similar effect is described as an higher-order effect 



m 



Busse However we stress here that this a-effect is in essence very different from the one that 



would be obtained with a flow geometry (i). 

Finally, we stress that the a-effect, so far understood as the generation of a toroidal mean electromotive 
force from a toroidal mean magnetic field, and described by only one a-coefficient, is in fact a special case 
of a more general connection between mean electromotive force and mean magnetic field, described by an 
a-tensor. We shall calculate the additional coefficients of this a-tensor and give their scaling properties 
in terms of the flow parameters like the number of rolls, the Rossby wave frequency and the magnetic 
Reynolds numbers (horizontal and vertical) of the flow. Depending on these parameters these additional 
coefficients may be dominant compared to the a-effect mentioned above and then completely change the 
overall picture of the possible dynamo mechanism. We shall see that it is all the more true for rolls not 
only drifting but also having a radial shift as expected for Rossby waves traveling in a spherical shell. In 
that case the rolls are bent as illustrated in figure [2)3. 

In section [2] we set the problem, give the basic equations, define the general assumptions and specify the 
velocity field. In section [3] we derive general analytical expressions for the mean electromotive force with 
respect to a cylindrical coordinate system. In section |4] we present asymptotic and numerical results for 
both cases (i) and (ii). Finally, in section [5] we discuss our results. 
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2 General concept 

2.1 Basic equations and assumptions 

We consider a body of a homogeneous electrically conducting incompressible fluid, penetrated by a mag- 
netic field B and showing internal motions with a velocity U. The magnetic field B is assumed to be 
governed by the induction equation 

f) V 2 B + Vx(UxB)-3 t B = 0, V • B = 0, (la,b) 
where rj is the magnetic diffusivity, considered as constant. 

Referring to a cylindrical coordi nate system (s, (p, z) we define mean fields as in Braginskii's theory of 
the nearly axisymmetric dynamo ( Braginskiil Il964a.bl ) by averaging over ip. Given a scalar field F, the 



corresponding mean field is denoted by F. It corresponds to the axisymmetric part of F. In the case of 
vectors or tensors the same definition applies to each component so that, e.g., B = (B s , B^, B z ). 
We split B and U into mean fields, B and U, and deviations b and u from them, that is, 

B = B + b, U = U + u. (2a,b) 

From the induction equation ([1]) we obtain 

??V 2 B |Vx(UxB) + Vx5-3,B = 0, V • B = 0, (3a,b) 

where 



e = u x b (4) 

is the mean electromotive force due to u and b. 

In view of the generation of a mean magnetic field, two terms of the equation ([3|) are of particular 
interest, that with U and that with £. We assume here that the mean velocity U corresponds to a rotation 
only. We further think of a proper specification of the small-scale velocity u so that £ covers the effect 
of Rossby waves. For the determination of £ we assume rigid-body mean rotation, that is we ignore any 
differential rotation, and adopt a co-rotating frame of reference in which U = 0. According to (|4]) £ is 
determined by u, which we consider as given, and b. Using (pQ), ([2]) and ([3]) we obtain 



r/V b + Vx(uxB) + Vx(uxb-uxb) — d t h = , V • b = . (5a,b) 

Clearly this equation determines b if u and B are given. 

In order to make analytical calculations possible we introduce a quasi-linear approximation (also known 
as the second-order correlation approximation in mean field theory), that is we neglect the term | V x (u x 
b — u x b) | in ([5]) . A sufficient condition for that approximation is 

min(iC S t ) « 1 (6) 

with the magnetic Reynolds number R' m = uol'/rj, and the Strouhal number St = uqt/1', and where no 
is a typical magnitude of u, /' a characteristic small length scale of the roll and r a characteristic time of 
the Rossby wave. If r is interpreted as the inverse wave frequency, St is the ratio of the roll's turn-over 
frequency to the wave frequency. At the end of section 12.21 we will come back on the condition ([6]) for the 
applicability of the quasi-linear approximation. 

We assume that the fluid velocity u is non-zero only inside a cylindrical layer with the mean radius Iq 
and the thickness 251$ (5 < 1), that is, in (1 — 8)Iq < s < (1 + 5)Iq. Its depen dence on s and (p and on time 



will be specified later. Moreover we consider u as independent of z (see e.g. iKim et Il999l for a similar 
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approximation). Clearly E can only be different from zero inside that cylindrical layer. As for B, which is 
by definition independent of (p, we assume that it is independent on z, too. Further its time dependence 
is considered as weak compared to that of u and therefore neglected in the following calculations. The 
independence of u and B on z suggests to consider also b as independent of z. Then also £ does not longer 
depend on z. Of course, in view of applications of our results to the Earth's core and the DTS experiment 
a treatment of the more general case with z-dependencies of u, B, etc. would be of high interest. This is 
left for future work. 

In what follows we measure all lengths in units of Iq, the time in units of Iq/tj and the velocity u in 
units of uq. For several purposes it is useful to split u into its parts and um perpendicular and parallel 
to the rotation axis. We measure and Uii in units of uq± and u u defined analogously to uq. With the 
assumptions introduced above, ([5|) turns then into 

V 2 b - d t b = - V x Q , V • b = (7a,b) 

with 

Q = R m uxB = (R m± u ± + R m ||U||) xB (8) 

and with the magnetic Reynolds numbers 

e> U l D D n 0||^0 , , 

K m — , %il — , K m \\ — . ^ya,D,cJ 

rj j] 11 rj 

In contrast to R' m , these quantities are defined with the large-scale parameter Iq. 
2.2 Fluid velocity 

Since the fluid is considered as incompressible, we have V • u = 0. As u is taken independent of z, we may 
represent it in the form 

u = — e x W + u z e , (10) 

where e is the unit vector in the z-direction. The stream function ip as well as u z may depend on s, ip and 
t. 

With the intention to simulate Rossby waves in their simplest form, we further specify the velocity u by 



ip = ip(s) cos(m(p — <p(s)) — ut), u z = u z (s) cos(m(p — j(s)) — ujt) (lla,b) 

with a positive integer m and constant u. The functions <p(s) and j(s) are the radial dependent phase 
shifts of the horizontal and vertical flow components. They can be enforced for example by the spherical 
geometry. In this c ase they correspond to rolls with horizontal section shapes like bananas (see for example 
Zhang et al. 20071 ) . Two examples of the flow streamlines are given in figure [2j Each pattern drifts in the 



azimuthal direction with the dimensionless angular velocity ui/m. 
For later purposes we define a complex vector u(s), depending on s only, such that 



u = Re{u(s) exp(— mup + iwt)} 



(12) 
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In a more explicit form, u(s) is given by 



an ~ 

u s = exp(imy?(s)), (13a) 

s 

u v = ~[d s ip(s) + imijj(s)d s (p(s)] exp(im<^(s)), (13b) 
u z = u z ( s ) exp(im7(s)). (13c) 



Various relations between the horizontal and the axial flow, that is between ip and u z , can be specified 
by the choice of <f(s) and j(s). As we are mainly interested in the two cases (i) and (ii) described in figure 
[U we define the two corresponding flow types, 

type (i) : ip(s) = j(s) , type (ii) : ip(s) = j(s) + vr/2m. (14a,b) 

The flow of type (i) corresponds to an axial flow driven by Ekman pumping. The zero lines of the axial 
flow coincide then with the borders of the cells of the horizontal flows. The flow of type (ii) corresponds 
to geometrical slope effect as the main cause of the axial flow. The extrema of the axial flow are then at 
the borders of the cells of the horizontal flow. 

Returning to the condition ([6]) for the applicability of the quasi-linear approximation we specify now I' 
and r such that I' = 51q and r = Iq/t/lj. Then we have 

R' m = uq51o/t] and St = uqIq/^Slo = R' m /5 2 uj. (15a,b) 

For stationary flows {to = 0), irrelevant for Rossby waves but still of interest for comparison with some 
numerical simulations, the condition (|6|), which is sufficient but not necessary for the validity of the 
quasilinear approximation, takes the form R' m <C 1. However, as turned out in simulations (Schrinner 
et al. 2005, 2006), this approximation may well apply for values of R' m up to the order of unity. For a 
non-stationary flow, if 5 2 lu S> 1 the condition ([6]) turns into St <S 1- In view of the Earth's fluid core, let us 
consider that Iq ~ 1800 km and r\ ~ lm 2 s _1 (which leads to l^/rj ~ 10 5 years). Assuming approximately 
equal extents of a roll in radial and azimuthal direction, 25 = ir/m, and m = 16 pairs of rolls we have 
5 ~ 0.1 (and then a radius 8lo of a roll about 180 km). A typical drift velocity urj/mlQ = 0.2 deg / year 
yields to u ~ 5.6 • 10 3 . Then we find that uj5 2 ^ 60. Then again a sufficient condition for the applicability 
of the quasilinear approximation reads St <C 1. It implies uq <C 0.3 mm s -1 , which gives the flow intensity 
upper-limit above which the quasi-linear approximation might not work. 



3 Calculation of S 

3.1 Poloidal-toroidal decomposition and reduction of equations 

In order to solve ([7]), we represent b as a sum of poloidal and toroidal parts, 

b = -V x (e x VS) - e x VT (16) 

with scalars S and T depending on s and cp. The components of b are then given by 



ldT dT 
b s = - — , b ip = — — , b z = -DS, (17a,b,c) 
s dip os 



where 



^ s ds^ S ds^ s 2 dip 2 ^ ^ 
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Likewise we represent V x Q in the form 

V x Q = V x (e x VF) + e x VG (19) 

with scalars F and G. Using ([7]), (|16p and (|19p we find, excluding singularities of F and G at s = and 

for s — > oo, 

DS-d t S = F , DT-d t T = G. (20a,b) 
With the help of the identity e • (V x (e x Vf)) = Df we further conclude from (fl~9j) that 

L>F = e-(VxQ), DG = e- (V x (V x Q)) = -e • (V 2 Q - V(V • Q)) . (21a,b) 
The first of these relations can be written in the form 



The second one is equivalent to DG = —DQ Z or, if we exclude again singularities at s = and for s — > oo, 
to 

G = -Q z . (23) 



3.2 Calculation of b 

Like the components of u, those of Q as well as the functions F and G have the form 

f(s, (f) = f^{s) cos(mtf - ujt) + f(*\s) sm(m(p - ut) = Re{/(s) exp(-im^ + iut)} , (24) 

where / is a complex quantity depending on s only. The same applies to b as well as S and T. In this 
notation the relations (|17p take the form 

b s = -—f(s), 6 =-^M, b z = -D m S(s) (25a,b,c) 
s as 

with 

The equations flU, ([22]) and ([23]) reduce to 

{D m -iu)S = F, (D m -iu)f = G, (27a,b) 



L» m # = -(|-(sQ V3 )+im4) J G = -Q z . (28a,b) 



In order to determine b we have to solve (I27p with i 7 and G satisfying (|28l) . Generic solutions of such 
equations are derived in Appendix According to (|A.8P - (|A.10p the solution of the first equation of (|28p 
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can be written in the form 

F(s) = -^^^[^(s%(s')) + imQ s (s')) s'ds', (29) 
where h m (s,s') is the Green's function defined in (|A.11|) . After integrating by parts, we have 

F(s) = J™ { dhm ^ S ' ] Q<p(s') - im^l^ Q s (s')) s> ds> . (30) 

Now from ([28]), ([H and with the help of (lATTT) - (lAT3ll . the solutions of ([27]) turn into 

n°° f)h - in? 

km(s,s'){^(s',s")Q v (s") - ^h^s'^^Qsis")) s"ds"s'ds\ (31a) 

/>oo 

f(s) = / k m {s,s')Q z (s')s'ds', (31b) 



J 

where k m (s,s') is the Green's function defined in (|A,4|) . 
Let us return to b as given by (|2"5j) . With (|3"Tj) and (f2"T|) we obtain 
f°° im , - 

b s (s) = - — A; m (s,s')Q^(s')s'ds', (32a) 
Jo s 

b< P (s) = -J ^W)Q*(*Vdy, (32b) 



, , , / / dh m (s, s ) a / /\ . h m (s, s ) a / /n\ / j / 
b z{s) = - I ( ^ (^(sj-im Q s {s))sds 



o 



oo />oo 



+ w / / fc m ( S ,/)(^(/, s , o^( S / o-^^( S ^ 5 '04( S / o)s ,, d S / vd S / . (32c) 



'0 JO 

After some algebra we further find 



MM / [( - 9h ™(s; s ^ +iuA 1 (s,s'))Q (p (s') + ^(h m (s,s') -iuA 2 (s,s'))Q s (s')] s' ds' (33) 





with 

roc o, / // / \ /-oo 

A 1 (s,s') = J k m (s,s") m{ d y >- s"ds", A 2 (s,s') = J k m (s, s")h m (s" , s')s"ds" . (34a,b) 

With the help of relations of Appendix [B] we can show that 

Id 1 

Ai(s,s') = — — (hm - k m )(s,s') , A 2 (s,s') = —(h m - k m )(s,s'). (35a,b) 

ILO OS \UJ 

This leads to 

k{s) =J°°[- 3km Q s ; S,) Q^s') + ^k m (s,s')Q s (s')] s'ds'. (36) 



May 5, 2009 10:12 Geophysical and Astrophysical Fluid Dynamics GAFD'Plunian'arxiv 

Rossby waves and ct-effect 

From the second relation ((7|) we have 



9 



Qs = R m ±u ip B z - R m \\u z B v , (37a) 
Qtp = Rm\\u z B s - R m ±_u s B z , (37b) 
Qz = R m ±{u s B [p - u v B s ) . (37c) 



Inserting this into (|32[) and ([36]) we obtain 



r°° k (<t <i') — — 

b s (s) = -imR m± / W ' > {u s {s') B v {s') - u^(s') B s (s')) s' ds', (38a) 
Jo s 

b v (s) = -R m± J ] {us(s')B^ S ')-u^s')B s (s'))s'ds', (38b) 

b z (s) = -R m \\ [°° [ dkm ^J S>) B s (s') + ^k m (s, 8')B v (s')] Ms') s' ds' 







S' 



+ R m ± J~ f^;^ u s (s') + '™k m {s,s')u 9 {s')}B z {s')s'ds'. (38c) 

As mentioned in Appendix[A] in the stationary case, that is to = 0, k m turns into h m . We note that as the 
fluid is at rest outside the moving layer, the integrations in (|38p and previous equations can be reduced to 
l-S<s<l + 6. 



3.3 Integral representation of S 
According to ([!]) we have 



yRe{u*xb}, (39) 



where u* is the complex conjugate of u. When inserting b s (s), b<p(s), b z (s) as given by ([38]) we see that £ 
can be written in the form 



£ K (s)= / K KX (s,s')B x (s')s'ds'. (40) 
Ji-s 



Here and in what follows k and A stand for s, ip or z. The kernel K K \ is then given by 



~ T Rm± \r J K * lf (KA) " {(ssUsvWsUw) ~~ (41) 
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The dimensionless quantities K K \, which do no longer depend on R m ± or R m \\, are given by 

AK ss (s, s') = -^r(*, s') u* v (s) u z (s') ~ ^(s, s')u* z (s) + ex., (42a) 

4K S(p (s, s') = — —k m (s, s') u* (s) Uz(s') + S ')K( S ) u s {s') + c.c., (42b) 

4K sz (s, s') = +^r(*, s') fi* (s) u 8 ( S ') + ™M*> s>) u* v (s) u^(s') + c.c., (42c) 

4K vs (s, s') = +— k m (s, s') u* z (s) u v (s') + -^r( s ; s ') w z (s') + c.c, (42d) 

~ 1772 

AK w (s,s') = — -k m (s,s')u* z (s)u s {s') + — k m (s, s') u*(s) u z (s') + c.c, (42e) 
4i^( S , S ') = ~^-(s, s') u* s ( S ) u s (s') - ™k m (s, s')u* s (s) u^s') + c.c, (42f) 

(I ft 1 TT) 

4K zs (s, s') = +~^r(s, s'K(s) u v {s') - —k m (s, s') u* v (s) u^(s') + c.c, (42g) 



4i^( S , s>) = ~^(s, s'K(s) u s (s') + ™k m (s, s') u* v {s) u s (s') + c.c, (42h) 



where c.c. denotes the complex conjugation. By symmetry reasons and because of the z-independence of 
the flow we always have K zz = 0. 

3.4 Expansion of S 

We rely now on the integral representation (|40p of £ K , assume that B\(s) varies only weakly with s and 
use the Taylor expansion 



In this way we obtain 



B x (s') = B x (s) + (s> - s) + ■■■. (43) 



£ K (s) = a KX (s) B x (s) + b KXs (s) + • • • (44) 

with 

rl+6 rl+5 
a K x(s)= K KX {s,s')s'ds', b KXs (s) = l K kX {s, s') (s' - s) s' As'. (45a,b) 

Jl-8 Jl-8 

The factors Iq have been inserted in (|44p and (|45p in order to give b K \ s the dimension of a magnetic 
diffusivity. The terms with higher derivatives of B\, indicated by • • • , are ignored in the following. 
Expressing K K \ in the relations (J45]) in terms of K K \ we obtain 

a K x = ^R m ± lo m 4 S «A if («A) = & ZS }4 Z<P ), w v , a KZ = 0, (46) 
Z l^mllJ V 7 \(ss),(sy?),(v?s),(^) 1 ; 



- f^W foj & k a, if («A) - { {ssUs(pU(psUtp(p) , " 0, (47) 
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where a K \ and b K \ s are dimensionless quantities independent of R m ± and R m \\ . Analytical expressions for 
the a K \ and b K \ s can be derived on the basis of (HT1) . (|42|) and ((131) . Of course, a 2 ^ = and b zzs = follow 
from K zz = 0. As for a sz = a vz = we note that K sz (s, s')s'ds' and f^f K ipz (s, s')s'ds' with K sz 
and K vz according to (142)) vanish. This can be shown with integrations by part and using V • u = 0. 

Let us add a remark on the nature of the expansion (|44fl and the coefficients a K \ and b K \ s . In most 
representations of mean-field electrodynamics the connection between £ , B and its derivatives is, with 
respect to a Cartesian coordinate system, given in the form Si = &ijBj + bij^dBj/dx^ + • • • . (Usually the 
notation a^- and b^ is used instead of and b^. We deviate from that, since a%j and b^ are already 
otherwise defined in this paper.) It is understood as a coordinate-independent connection, which implies 
that Si and Bj are components of vectors, and Sjj, as well as dBj/dx^ components of tensors, all 
with the well-known behavior of such objects under coordinate transformations. The a K \, however, do 
not completely coincide with the components of the tensor derived in that sense from the a«. The reason 
is that the transformation of dBj/dxk in our cylindrical coordinate system produces not only terms 
with derivatives of B K but also such without derivatives (the same remark would apply if our coordinate 
system was spherical instead of being cylindrical). In the common understanding the a-effect is, again in 
coordinate-independent manner, defined on the basis of the contribution dijBj to Si. We slightly deviate 
from this definition in what follows. When speaking of a-effect we refer simply to the contribution a K \B\ 
to S K . This is in so far justified as it is just this contribution which describes, e.g., the generation of S v 
from Bio. 



4 Some typical examples 

4.1 Specification of the flow patterns 

In this section we present numerical results for some typical flow profiles corresponding to rolls of type (i) 
and (ii) as defined in (jl4() . In addition we also consider two types of phase-shift radial dependence, 

(a) 7(a) = and (6) j(s) = ^ + / ~ ^ . (48a,b) 

o 

The flow geometry is further specified by 

n, = ^(i-a 2 , t = 5(i-ef, e = ^, if iei<i, 

u z = ip = otherwise . (49) 

Here u z and ip are normalized such that at any time the average of u z over a surface given by 1—5 < s < 1+5 
and — 7r/2m < ip < it /2m as well as the average of at a given value of (p over 1 < s < 1 + 5 are equal 
to unity. Figure [2] shows isolines of u z for the two case (a) and (b) defined in f|48j) . 

In case (a) the extent of a flow cell is 25 (in units of Iq) in the radial direction and ir/m in the azimuthal 
direction. We speak of "compact rolls" if their ratio 25m/ ir, or simply 5m, is in the order of unity. In that 
sense Fig. (2^i shows compact rolls. In case (b) the radial phase-shift leads to extended rolls as shown in 
Fig. [2b even if 5m is in the order of unity. For simplicity we always set 25m/'K = 1 in the rest of the paper. 



4.2 A first analysis of the results 

The results obtained so far allow us to draw some conclusions concerning the structure of E and the 
a-effect. Let us consider the kernel K for flows of types (i) and (ii) defined in (j!4p and recall that their 
axial parts are driven by Ekman pumping or the geometrical slope effect, respectively. 
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In the stationary case (u = 0), irrelevant for Rossby waves but still of general interest, with the 
help of (|14p and (|42p . we conclude that K has the following structure 



K = X X ] in case (i) , K = X x ) in case (ii) , (50a,b) 





where crosses stand for matrix elements which are not necessarily equal to zero. In case (i) we see from 
(|50D that K lea ds to non-zero coefficients a ss and a w corresponding to the dynamo mechani sm described 



bvTBussel (I1975T ) which can be interpreted as an a -mechanism. In the Karlsruhe experiment (jMuiler et all 



20041 . 120061 7 the dynamo action of a mechanism of that kind has been demonstrated. On the other hand 



in case (ii) ([50]) implies that a w = (and a ss = 0) in accordance with the heuristic arguments of figure [TJ 

For drifting waves (uj ^ 0), the kernel components except K zz are not necessary equal to zero. 
Then it may happen that the coefficients off-diagonal be dominant and therefore the simple a-effect 
mentioned earlier not be a relevant part of the dynamo mechanism. In fact, in addition to the a-effect a 
transport of mean magnetic flux by the so-called 7-effect can be expected. To derive it we can define a 
symmetric matrix a and a vector 7 in the following way 

a K \ = -^(a K A + a\ K ) , 7 K = 2 e «V a V> (51a,b) 

where e K \^ has to be interpreted in the usual way (Levi-Civita symbol) identifying the subscripts s, (p and 
z with 1, 2 and 3, respectively. Then we can write the mean electromotive force in the form 

£ = -aoB-7xB (52) 

in which the s-derivatives of B are not included. As can be seen from the structures of the matrix [a K \] 
contributions to both a and 7 are expected. If the role of the coefficient a w (usual a-effect) in the dynamo 
process is well understood as being directly related to the generation of a poloidal field from a toroidal 
field, the role of the other coefficients of a is less clear and would need a specific study in itself. 



4.3 Numerical results 

In this section we plot the coefficients a K \ for flows of type (ia) in figure [3l (iia) in figure [H (ib) in figure 
[5] and (iib) in figure [6l We consider rolls satisfying 8m = tt/2 and show results in the double asymptotic 
limit 5 <^ 1 and 5 2 uj 1. We find that the profile of each coefficients a K \ converge in this double limit. 
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We also find some scaling laws in 5 and u such that 

~a w ~ uPS". (53) 

We found that their validity prevails even for uj5 2 ~ 1. For flows of type (ia) and (iia) the scalings in oj is 
consistent with the general structure of the kernel K given in (|50|) for uj = 0. 

For flows of type (ib) and (iib) the scalings are rather different from cases (ia) and (iia). In particular we 
found no uj dependency, suggesting that the leading term in k m (s, s') is of order 8. Surprisingly enough in 
case (iib) we found that a w ~ ruling out any chance to explain the dynamo mechanism with a simple 
a-effect. 

It is interes ting to compare the results obtained for the a-effect coefficient aj_ in the case of the 
Roberts flow (jBadler et all l2002ah to our results for a w for a flow of type (i) with compact cells. 



The length a used there corresponds to 25lo, and the Reynolds numbers Rm± and Rm\\ used there 
have to be interpreted as 25R m ± and 5R m \\, respectively, with our R m ± and R m \\. In that sense the 
result a± = (ir 2 /16)(rj/a)Rm±Rm^(p(Rm±) reported in the mentioned paper takes the form aj_ = 
(it 2 /16)(r]S/lQ)R m ±R rn \\(j)(R m ±). This quoted result applies to arbitrary Rm± and Rm\\. The function cj) is 
equal to unity for Rm± = and so in the second-order correlation approximation. It decreases monotoni- 
cally if Rm± grows and tends to zero as Rm±_ — > oo. According to (j46j) we have a w = {r]/lo)Rrn^R m \\0'w 
The same relation applies with the averages < a w > and < a w > of a w and a w over s. We take from 
figure [3] that < a w >ss 5. This leads to < a w >ss {tl^ /h)Rm±Rm\\ ■ Hence our result for a w (derived 
in the second-order correlation approximation) is in reasonable agreement with the result for the Roberts 
flow. 

The comparison with the result for the Roberts flow suggest that our result for a w remains valid for all 
values of R m \\ and, as can be concluded from the specific properties of 4>, for values of <5i? m || up to the 
order of unity. There is, however, no straightforward extension of the proof for the linearity in i? m || to our 
case. We further learn here that it is of less importance for the magnitude of the a-effect whether or not 
a given roll is at all sides surrounded by other rolls. 



5 Conclusions 

The main insight of this study is that the a-effect and presumably the dynamo action generated by Rossby 
waves depends drastically on the Ekman number E. The reason for that lies in the competition between 
Ekman pumping and geometrical slope effect as the main source of the axial flow, producing different phase 
shifts between the horizonthal and vertical components of the flow. In the limit E <C 1 the geometrical 
slope effect prevails. In particular we find that the mean-field coefficients have then different radial profiles 
and, probably more important, different asymptotic scalings in the limit of small rolls radius 4 < 1 and 
wave frequency such that 5 2 oj <C 1. 

In addition we found that for rolls with an arbitrary radial shift the a-tensor is completely changed not 
only in terms of the radial profile of the tensor coefficients but also in terms of scalings in 5 and to. In 
particular we found that the a-effect corresponding to the coefficient a w may disappear, making then 
difficult the interpretation of the dynamo mechanism if any. 

Th ere are two recent numerical studies in which the a-tensor produced by Rossby waves has been calcu- 
lated. ISchrinner etdl (1200,1 [2006) extracted the a-tensor (and higher-order coefficients) from geodynamo 



simulations with E = 10 3 . Schaeffer and Cardin ( 20061 ) did the same with simulations of the Taylor 



Couette convection using the quasi-geostrophic approximation with values of E down to 10~ 8 . It is then 
of interest to compare qualitatively our results to these findings (ignoring the difference addressed at the 
end of section [3741 between the two defini t ions o f the a-effect a w and a w ). 
In the study by ISchrinner et al. I (|2005l . l200fil ) a thermally driven dynamo in a rotating spherical shell 



is considered. In a regime not too far from to the onset of convection, drifting rolls are observed. As a 
consequence of E = 10" 3 , the main cause of the axial flows in these rolls should be Ekman pumping. 
Indeed the dependence of a w on s corresponds roughly to that of our d w depicted in Figures [3l This 
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Figure 3. Scaling and s-profile of coefficients a K \, for 5m = 7r/2, <5 1 and <5 2 o; <C 1, and for rolls of type (ia). 
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Figure 4. Same caption as figure [3] but for rolls of type (iia). 
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Figure 5. Same caption as figure [3] but for rolls of type (ib). 
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shows that for not a too small Ekman number, the dynamo process can be understood on the basis of the 
mechanism described in figure [TJi) . 
Schaeffer and Cardinl (2006) dealt with the Rossby wave instabilities of a geostrophic internal shear layer 



produced by differential rotation between two spheres. They used the quasi-geostrophic approximation 
including both Ekman pumping and slope effect as sources of the axial flow. By inserting the corresponding 
velocity field in an induction equation solver they showed that it is capable of dynamo action. They also 
derived the a coefficients from their numerical simulation for E = 1CP 8 . The dependence of their a w on 
s corresponds roughly to that of our o w when adapting our flow definition (|49j) to the quasi-geostrophic 
approximation. The comparison of the other coefficients is left for future work. 
Let us add a remark concerning the sign changes of a w in the case (iia). Investi gating a spherical a 



effect dynamo model with isotropic a-effect and a spherically symmetric coefficient a. lStefani and Gerbeth 
(200 5|) found polarity reversals of the mean magnetic field if a changes its sign along the radius (see also 
Giesecke et Ao05al lbh. 



In our calculations both the fluid flow and the magnetic field have been considered as independent of 
z. Therefore the results reflect by far not all essential features of the convection rolls in a rotating liquid 
sphere. A more realistic treatment of the problem would require to consider not only the z-dependencies 
mentioned but also the multi-scales structure of the flow. The difficulties that arise in this way could 
perhaps be reduced by taking advantage of the quasi-geostrophic approximation. It is then possible that 
only rolls of a certain scale are important for the dynamo process, that is, those having a sufficiently 
large 5R m to produce dynamo action but not too large in order to avoid too strong flux expulsion and 
a-quenching. In this case the simple picture of a ring of drifting rolls generating an a-effect as one part of 
the dynamo mechanism (completed by differential rotation) might be still relevant. 
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Appendix A 

General solutions of equations (|27p and (|28p in terms of Green's functions 
In view of (|27p we consider first the equation 

(D m - iuj)f(s) = g(s) (A.l) 

in < s < oo. Its general solution reads 

f(s) = c 1 I m (Vkj s) + c 2 K m (yAuJs) 

— K m (\f\bj s) / I m (V\w s') g(s', k) s' ds (A. 2) 

Jo 

/oo 
K m (Viuj s') g(s', k) s' ds' , 

where c\ and C2 are arbitrary constants, and I m and K m modified Bessel functions of the first and second 
kind (see Appendix [B]) . In order to exclude singularities at s = and for s — ► oo, we use (|B.2p and (|B.3p . 
put c\ = C2 = and obtain 



/•oo 

f(s) = - k m (s,s')g(s')s'ds', 
Jo 



(A.3) 
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with a Green's function k m denned by 

k m (s, s') = K m (Vhv s)I m (v / iw s') for s' < s 

k m (s, s) = I m (V\uj s)K m (Viuj s') for s < s' . (A. 4) 

Applying this to P7|) . we find 

poo poo 

f(s)= k m (s,s')Q z (s')s'ds', S(s) = - k m (s,s')F(s')s'ds'. (A.5) 
Jo Jo 

The Green's function k m defined by (|A.4p is complex. Its real and imaginary parts and can be 
written in the form: 

k^(s,s f ) = ker m (v / LJs) ber m (-\/~ujs') — kei m {-\fujs) bei m (-\fujs') for s' < s 

= ber m (\fujs) ker m (v / ws / ) — bei m (y/ujs) kei m (-\/ujs') for s' > s (A. 6) 

A;^(s,s') = ke.i m {\fu)s) ber m {^/u)s'^ +ker m (v / ws) bei m {yfujs') for s' < s 

= bei m (\/iHs) ker m (y / ws / ) + ber m (y/ujs) kei m (-\/ujs') for s' > s (A. 7) 

where ber m , bei m , ker m and &ei m are the Kelvin functions of order m. 

In view of (|28[) we consider first the general ordinary differential equation 

D m f(s) = g(s) (A.8) 

in < s < oo. Its general solution reads 

/(,) = Cl s m + c 2 s~ m - ±- (s~ m J g(s') s' m+1 ds' + s m j g(s') s'~ m+1 ds') , (A.9) 

where c\ and C2 are again arbitrary constants. In order to exclude singularities at s = and for s — > oo 
we put ci = C2 = 0. In this way we obtain 

/•oo 

f(s) = ~ h m (s,s')g(s')s'ds' (A.10) 
Jo 

with a Green's function h m 

hm(s,s') = ^(^rfor s ' < s 

Ms,s') = ^(4r fOTS < s/ - (A. 11) 

As to be expected, k m turns into h m if u; = 0. 

Appendix B 

Useful relations for modified Bessel functions 

The modified Bessel functions of first kind, I m (x), and of second kind, K m (x), are solutions of the equation 

x 2 y" + xy — (x 2 + m 2 )y = 0, m > 0. 
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From the recurrence relations 

in 

l' m {x) = +I m -i{x) - —I m (x) 

Tfl 

K L( X ) = -K m -i{x) - —K m (x) 

and the Wronskian relation 



I m {x)K rn „i(x) + I rn ^i(x)K m (x) = I m {x)K m+ i(x) + I m+1 (x)K m (x) = ^ 



we conclude that 



I m (x)K' m (x) - l' m (x)K m (x) = ~ . (B.l) 



The modified Bessel functions satisfy the integral relations 

J x m+1 I m (x)dx = x m+1 / m+1 (x) , J x- m+1 I m (x)dx = x- m+1 7 m _i(a;) 
J x m+1 K m (x)dx = -x m+1 K rn+l {x) , J x- m+1 K m {x)Ax = -x' m+1 K rn ^(x) 



as well as the asymptotic relations 



and 



l m (x)«i r (fr(l + ^(f)») as x^O 

Kl (x)^±(l + ^ log(|)) as x^O (B.2) 



i nr 

I m (x) ~ e x , K m (x) m \ — e x as x — > oo. (B.3) 

V27TX V 2x 



